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Abstract 



A model of aluminium has been developed and implemented in an Ocean 
General Circulation Model (NEMO-PISCES). In the model, aluminium enters 
the ocean by means of dust deposition. The internal oceanic processes are 
described by advection, mixing and reversible scavenging. The model has 
been evaluated against a number of selected high-quality datasets covering 
much of the world ocean, especially those from the West Atlantic Geotraces 
cruises of 2010 and 2011. Generally, the model results are in fair agreement 
with the observations. However, the model does not describe well the vertical 
distribution of dissolved Al in the North Atlantic Ocean. The model may 
require changes in the physical forcing and the vertical dependency of the 
sinking velocity of biogenic silica to account for other discrepancies. To explore 
the model behaviour, sensitivity experiments have been performed, in which 
we changed the key parameters of the scavenging process as well as the input 
of aluminium into the ocean. This resulted in a better understanding of 
aluminium in the ocean, and it is now clear which parameter has what effect 
on the dissolved aluminium distribution and which processes might be missing 
in the model, among which boundary scavenging and biological incorporation 
of aluminium into diatoms. 
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The distribution and cycling of aluminium (Al) in the ocean has received atten- 
tion for a variety of reasons. Firstly, if the Al cycle is understood well, aluminium 
surface concentrations can be used to constrain atmospheric dust deposition fields 
(I Gehlen et ah] (|2Q03[ ); [Measures et~al] ( |2QQ5[ ); |Han et al.| (|2QQ8|); [Measures et~ 
(2010); Han ( 2010[ )), which are used to predict aeolian iron addition to the euphotic 
zone. This is important, since iron is an essential trace- nutrient for phytoplankton, 
thus its availability has a direct consequence on primary production and air-sea 
CO2 exchange. Secondly, there is evidence that Al inhibits the solubility of sedi- 
mentary biogenic silica (jLewin (1961); Bennekom et al. ( 1991| ); Dixit et al. (2001); 



Sarmiento and Gruber (2006) 



Emerson and Hedges (2006)). If less biogenic silica 



gets dissolved from sediments, eventually there will be less silicic acid available in 
the euphotic zone, which will reduce diatom production as silicon is an essential 
major nutrient for diatoms. Modified diatom productivity will impact ocean food 
webs and the export of organic carbon to the ocean's interior. For advancement in 
both of these fields of interest a good understanding of the Al cycle is pertinent. 
Currently it is assumed that the major source of Al to the ocean is via dust 



deposit ion (e.g. jOrians and Bruland[ ( [1986[ ) ; [Maring and Duce[ ( |1987[ ) ; [Kramer et al. 
( 2004| ); Measures et al. (2005 )). When dust enters the ocean, a part of its aluminium 
content (1-15%) dissolves in the uppermost layer and is quickly distributed over the 
mixed layer by turbulent mixing. Most dust remains in the particulate phase and 
sinks to the bottom of the ocean, while a small fraction might dissolve in the water 
£olumn. The dust that does not dissolve at all is buried in the sediment ([Gehlen 
|et al.[ ( [2003| ); [Sarmiento and Gruber] ( [2006| )) 



Arguments that dissolution occurs primarily in the upper layer of the ocean come 



from shipboard experiments and atmospheric moisture considerations. IMaring and 



puce 



1987) and Measures et al. (2010) showed that within a day after deposition, 
most of dissolvable Al will be dissolved. Assuming a sinking speed of dust of 30 
m/day, most Al would then dissolve in the upper 30 m of the ocean. This depth is 
shallower than the mixed layer depth, which means that there is little dissolution 
below the mixed layer. Even though some earlier studies showed that most deposi- 
tion is dry (Jickells et al. (1994); Jickells ( 1995|)), more recent work shows that dust 
deposition is mostly wet (Guerzoni et al. ( 1997| ); [Vink and Measures ( 2001| )). It 
has been argued that most dust passes through low pH environments in the atmo- 
sphere, which means that for wet deposition Al is already dissolved when it enters 
the ocean surface. Since the wetly deposited Al is most important, dissolution in 
the surface ocean is most relevant (as a lower bound, since dry deposition results 
in both surface and water column dissolution). Furthermore, based on the same 
low pH argument, dissolvable Al from dry deposition is likely to instantaneously 



dissolve in the surface ocean (Measures et al. ( 2010) and references herein), making 
the relative amount of surface dissolution even higher compared to water column 
dissolution. 

Fluvial input can be thought to be important as well, since rivers carry large 
concentrations of Al, but in estuaries and coastal regions this Al is removed by 
scav enging of Al onto particles ( Mackin| ( [1986 ); Mackin and Aller (1986); Orians 
pmd Bruland] ( [1986[ ); [Brown et al.| ( |2010[ )). There are also indications for Al in- 



jputas a consequence of sediment remobilisation, as in the Arctic Ocean (Middag 
|et al.[ ( [2009l )) and North Atlantic Ocean ( jMoran and Moore[ ( |T99T| ) . However, the 
importance of sedimentary sources can vary by basin (e.g. in the Southern Ocean 
these are small as shown by Moran et al. (1992) and Middag et al. (2011)). Finally, 
hydrothermal vents are thought to only play a minor role (e.g. Hydes et al. (1986); 
Middag[ ( [2QTQl )). In summary, the dominant external source of aluminium in the 



ocean is atmospheric dust deposition. 
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p Dissolved aluminium (Aldiss) is removed mainly by particle scavenging ( [Orians 
[and Bruland] ( |1986[ ); poran and Moore| ( |1989[ ); [Bruland and Lohan] ( [20061 )). This 
is the combination of adsorption onto a solid surface, followed by sinking due to 
^Insufficient bu oyancy of the particulates in the seawater ([Goldberg (1954); [Bacon 
[and A nderson ( 1982 ) ; pruland and Lohan ( 2006| )). Typically, scavenging is deemed 
to be reversible, which means that during sinking release of the adsorbed, or par- 
ticulate, aluminium (Alads) occur. This happens both directly (by desorption) 
and indirectly (by dissolution of the biogenic carrier particles). As a consequence, 
Aldiss concentrations increase with depth (Bacon and Anderson (1982J ; Anderson] 
( [2QQ6| ). In this way Aldiss is distributed over depth more efficiently than due to 
mixing and water mass transport. Aluminium is scavenged relatively efficiently and 



therefore has a relatively short residence time in the ocean (100-200 yr) (Orians 



and Bruland (1985)). 



Except for scavenging, there are strong suggestions from observations in certain 
regions that Al is biologically incorporated into the siliceous cell walls of diatoms 
dStoffyn] ( [19791 ); |Moran and Moore] ( |1988D ; [Gehlen et al.| ( |2002D ). It seems that 
Al does not play an essential role for the diatoms, but it can be incorporated 
functioning as a replacement for silicon (Si), since it is similar in size. Therefore 
it is likely that the incorporation ratio Al : Si is close to that of the surrounding 
waters. These regions include the Arctic Ocean (Middag et al. (2009)) and the 



Mediterranean Sea (MacKenzie et al. (1978); Hydes et al. (1988); Chou and Wollast 
( 1997| )). Given the ratio Al : Si of incorporation into the diatom in the photic zone, 
after remineralisation anywhere in the water column, the same dissolved Al : Si will 
be present, as long as this is the only source of Al and Si. When the dissolved Al 
and Si is then advected into the Atlantic Ocean by the North Atlantic Deep Water 
(NADW), this signal slowly disappears because of other sources of Al among which 



dust deposition and possibly sediment resuspension (Middag et al. (2011)) and a 
source of Si from Antarctic Bottom Water. 

Recent years have seen the development of models of the marine biogeochemical 
cycle of Al 



Gehlen et al. (2003) implemented a basic scavenging model, while 



Han 



^t al. (2008) also included a biological aluminium incorporation module. Gehlen 
al. (20031) had the objective of testing the sensitivity of modelled Al fields to 
dust input and thus to evaluate the possibility for constraining dust deposition via 
Aldiss- To this purpose they embedded an Al cycle in the HAM0CC2 biogeochemical 
model. The model consists of an equilibrium relation between, on the one hand, 
Alads and, on the other hand, Aldiss- In chemical equilibrium Alads is proportional 
to the biogenic silica (bSi02) concentration. In their work, as well as this paper, 
the term biogenic silica or bSi02 refers to the detrital fractions which is fuelled by 
diatoms and other silicifying phytoplankton, which have no stable organic matter 
coating and sink. When bSi02 sinks to the seafioor (together with adsorbed Al), it is 
buried. The resulting concentration of modelled Aldiss was of the same order as the 
then published observations, but it suggested a significant overestimation of Saharan 



dust input (Gehlen et al. (2003)) when the dust deposition field of Mahowald et al. 
( [19991 ) was used. The main goal of Han et al.M2008k was to better constrain the dust 
deposition field. For this purpose they used the Biogeochemical Elemental Cycling 



(BEG) model improved by Moore et al. (2008) as a starting point. They used all 



dissolved Al datasets used by Gehlen et al. (2003) and added more datasets. Except 



for scavenging Han et al. (2008) added a biological Al uptake module where the 



Al : Si uptake ratio is a function of the ambient Al and Si concentrations (Han et al. 



( 2008| )). The residence time of Al for both modelling studies varies strongly between 



different locations (from less than one year to almost 200 years), but generally the 
^values ar e significantly lower than early estimates of 100-200 years (Orians and 
|Bruland| ([T985)). 

Overall, there are a number of questions regarding the oceanic Al cycle that re- 
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main to be fully addressed. These touch on issues of ocean circulation, the specific 
sources and sinks of Al in different parts of the world ocean and what processes 
are needed to accurately simulate the oceanic distribution of Al. Firstly, there 
is the question of the meridional (north to south) distribution of Al through the 
Atlantic Ocean. In the North Atlantic Ocean and northern seas, water sinks and 
forms NADW, which is then transported southward (e.g., Lozier ( 2010) ; Gary et al. 
(2011)). In the deep Atlantic Ocean dissolved silicon (Si) concentration increases 



from north to south (Ragueneau et al. ( 2000 ) ; iSarmiento and Gruber ( 2006[ )), while 
the concentration of Aldiss stays relatively constant until about 20° S and then de- 
creases (R. Middag, in preparation; see also Section [3|. Thus it has a generally 
opposite behaviour compared to Si. Since there are strong suggestions that the 
processes controlling the distribution of Si and Al are linked, the question is raised 
how this negative correlation is possible. Secondly, there is the question about the 
observed profiles of Aldiss at different locations in the ocean. Generally, profiles of 
Aldiss have a reversible-scavenging profile (increasing with depth) and often with a 
minimum near 1 km depth and a maximum at the surface because of dust deposi- 



ion. Howev er, observations in the Mediterranean Sea (Hydes et al. (1988); Chou 
■pind Wollast|(|1997) ) and IPY-Geotraces-NL observations in the eastern Arctic (Mid- 



dag et al. ( |2009[ )) show that there is a strong positive relation between aluminium 
and silicon. This supports the hypothesis of biological incorporation of aluminium 



into the ceh wah of diatoms ( [MacKenzie et al.| ( [1978) ; Stoffyn] ( |1979| )). 

These issues can be analysed further by the use of numerical models. Since there 
is a strong spatial variation in aluminium concentration (and its relation to silicon), 
an ocean general circulation model should be used to simulate the distribution of 
Aldiss- Potentially crucial parameters and sources can be modified in the model 
to test its sensitivity to these changes. In this way a better understanding of the 
aluminium cycle can be reached. 

In this paper we try to model the observed distribution of Al and examine the 
processes driving it. Based on new observations and previous work on aluminium 

( [20081 )) we formulate a model of alu- 



modelling (Gehlen et al. (2003) 



Han et al. 



minium based on dust deposition and scavenging by biogenic silica. This model and 
the configuration of the simulations will be set out in the following section. Then 
the observations which we used to check and improve the model will be discussed. 
The results of the several experiments follow in Section [Sj as well as a comparison 
with the observations. The discussion in Section [4] comprises of a comparison be- 
tween our model results and Gehlen et al. ( 2003| ), a timescale analysis and based on 
^hai a discussion of our simulations. We have not compared our results with |Han 
et al. (2008), since we have not performed simulations with biological aluminium 



incorporation. The last section ^ sums up the conclusions, gives an outline for 
further development of the model and suggests what further study is needed. 



2 Methods 
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2.1 Model description 

To model the three-dimensional distribution of Aldiss ^ the Ocean General Circulation 
Mode l (OGCM) cal led Nucleus for European Modelhng of the Ocean (NEMO) is 
used ( |Madec| ( 2008[ )). For this study we use PISCES, one of the biogeochemical 
components available in NEMO ( |Aumont and Bopp| ( |2006| ); |Ethe et al.] ( |2006[ )), 
which has been employed for many other studies concerning trace metals, as well as 
large scale ocean biogeochemistry (Aumont and Bopp (2006 ); |Dutay et aL (2009); 
Arsouze et aT] ( |2009[ ); |Dutay et al] ( |2009[ ); [Tagliabue et al.| ( |2010[ )). PISCES is run 



off-line forced by a climatological year of physical fields including subgrid turbulence. 
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20 60 100 140 180 



Figure 1: Velocity at surface, yearly average, on the 0RCA2 model grid. The colour 
represents the speed in m/s, while the vectors represent the direction of the flow. 
Vectors are plotted at every fourth gridpoint. 
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The forcing frequency was set to five days. These physical fields were calculated 



by the physical component of NEMO named OPA (Ma dec et al.| (|1998|), which 
was forced by satellite derived wind stress data. Details on the forcings and model 



configuration can be found in Arsouze et al. (2009) and Dutay et al. (2009). 



All input and output fields are defined on the 0RCA2 grid, an irregular grid 
with a nominal resolution of 2° x 2°, with the meridional resolution increased near 
the equator and two 'north poles' in Canada and Russia to eliminate the coordinate 
singularity in the Arctic Ocean. In Fig. [l] the annually averaged surface velocity 
field of the OPA output is plotted on the model grid. It shows the main features of 
the ocean circulation such as the equatorial current systems, the western boundary 
currents (Gulf Stream and Kuroshio) as well as the Antarctic Circumpolar Current. 



2.1.1 Biogeochemical model 
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We implemented the aluminium cycle (see Section 2.1.2) to PISCES, a biogeo- 
chemical model simulating the cycle of carbon and the major nutrients (nitrate, 
phosphate, ammonium, silicic acid, and iron), along with two phytoplankton types 
(diatoms and nanophytoplankton) , two zooplankton grazers (micro- and meso- 
zooplankton), two classes of particulate organic carbon (small and large) of dif- 
ferential labilities and sinking speeds, as well as calcite and biogenic silica. The 
standard version of PISCES accounts for 24 tracers. For a more detailed descrip- 



tion of PISCES see the auxiliary material of Aumont and Bopp (2006). 



The aluminium model interacts with the silicon cycle. The model PISCES dis- 
tinguishes three silicon pools: the silicon content of living diatoms (diatom Si); the 
silicon content of dead, sinking diatoms (biogenic Si) and dissolved silicic acid. In 
the model silicic acid and other nutrients are supplied to the ocean by means of 



atmospheric dust deposition, rivers and sediment mobilisation (Aumont and Bopp 



(2006); Gehlen et al. (2007)) 



185 
186 



2.1.2 Aluminium model 

For the modelling of Al, we introduced two new tracers to PISCES: dissolved alu- 
minium (Aldiss) and adsorbed aluminium (Alads)- Our Al model follows the ap- 
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proach of Gehlen et al. (2003), i.e., the internal processes (described below) of ad- 
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sorption and desorption are identical to the model of Gehlen et al. (2003), although 



the associated biogeochemical model is more complex. 

For the standard model configuration, the sole source of Al to the ocean is by 
dissolution of dust particles in the ocean surface. In one of the sensitivity exper- 
iments, a small fraction is also dissolved in the water column (see Section 2.1.4[ ). 
We assume an Al : Fe dust fraction of 8.1 : 3.5, which is based on the mass percent- 
ages of Al and Fe known to be present in the Earth's crust (Wedepohl ('1995')). 
The solubility of Al from dust is not well constrained, with reported values ranging 
from 0.5-86%, but is probably in the range of 1-15% (Orians and Bruland ( 1986[ ); 
Measures and Vink (2000); Jickells et al. ( 200 5*)). Furthermore, it is likely that th e 
source of dust has an impact on the solubility of the Al fraction (Baker et al. (2006); 



Measures et al. (2010) and references herein). We take the dissolution of dust as 



a constant with respect to the origin and location of deposition. The dissolution 
occurs only in the upper model layer, and is described by the following equation: 



d_ 
di 



[Al, 



■dissj 



(1) 



deposition 



where Vai is the Al flux into the ocean, a is the fraction of Al that is dissolved, and 
Hi is the thickness of first model layer. Square brackets denote the concentration 
of the tracer. Since in the mixed layer has rapid vertical mixing (z^^; is very large 
in the ocean model; the typical mixing time is Dml/^e)^ adding Aldiss in the first 
layer is effectively the same as adding it spread out over the mixed layer. The Al 
that does not dissolve from dust, either at the surface or in the water column, is 
assumed not to play any role in the biogeochemical cycle of Al and on our timescales 
of interest, and can be thought of as being buried in marine sediments. 

Dissolved Al is assumed to be of the form of aluminium that may adsorb onto 
biogenic silica particles and, aside from external inputs, the Aldiss concentration is 
essentially governed by adsorption and desorption. The Aldiss and Alads concentra- 
tions follow the following reversible first-order adsorption equation (Gehlen et al. 
( |2003J ): 



d_ 
dt 



[Al, 



ads J 



^([Al 



ad/desorpti( 



adsJ 



[Al. 



■adsJ 



(2) 



where 



IKl] = ■ [Aldiss] • [bSiOa] 



(3) 

in which bSi02 stands for biogenic silica. [Al^^^] is the chemical equilibrium con- 
centration of Alads- The parameter kd is the partition coefficient and n is the first 
order rate constant for equilibration of Alads to Al^^^. Since total Al is conserved 
when only internal processes are concerned, the time derivative of Aldiss is equal to 
the negative of the time derivative of Alads- 

Eq. [3] describes the chemical equilibrium between Aldiss and Alads 7 and Eq. [2] il- 
lustrates that some time is needed before equilibrium is reached (it is not modelled 
as an instantaneous process). As noted by Gehlen et al. ( 2003| ), this is a purely 
empirical model and a more mechanistic description of the adsorption/desorption 
of Aldiss onto bSi02 is difficult, since in reality, the pool of observed "dissolved 
aluminium" is operationally defined, i.e., it exists as different kinds of ions as well 
as colloids that can go through a filter of, e.g., 0.2 jiia. Also the adsorbent, bSi02, 
exists in various forms and it is, for instance, difficult to even define the charge 
density, which is important for its potential to adsorb Aldiss- Usually the surface 
is covered by natural organic matter (Lead and Wilkinson ( |2007 ); Filella ( 2007[ )). 
Furthermore, the main scavengers for aluminium in the open ocean also include Par- 
ticulate Organic Matter (POM), calcium carbonate (CaCOs) and even particulates 
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that are not part of PISCES (such as hthogenic particles that do not participate 
actively in biogeochemical cycles). Hence it is difficult to conceive a mechanistic 
model of ad- and desorption of "Aldiss" onto the poorly-defined pool of particulates. 

In our model we use only bSi02 as scavenger of Aldiss- Biogenic silica is likely to 
be the most important because (1) it sinks quickly, so that Al has little chance to sig- 
nificantly desorb from bSi02 (ballast effect), and (2) positively charged aluminium 



ions are easily adsorbed onto the negatively charged surface of bSi02 particles (Lou- 



Icaides et al. 



f2010')). While particulates are also presumed to be negatively charged, 
observational studies confirm that biogenic silica is a strong scavenger of aluminium 
dStoffyn and Mackenzie] ( | 1982 1 ); |Mor an and Moorej ( |1989| )). 

Sinking of bSi02 and Alads is described as follows. In the model, bSi02 sinks 
with a speed of wml = 30 m/day in the mixed layer. Below the mixed layer, it 
increases linearly with depth, reaching a value of w^]^^ = 200 m/day at 4 km below 



the mixed layer (Aumont and Bopp (2006')). This means that Alads also sinks with 



this speed, while it can desorb following equations [2] and [3] For every layer in the 
model, we have the following equation: 



d_ 

m 



[Al 



■adsj 



sinking 



d_ 
dz 



{Ws ■ [Alads]) , 



(4) 



where Ws is the sinking speed, given in the mixed layer by a constant wml and 
below the mixed layer increasing with depth according to: 



Ws ^ML + (^4km - ^Ml) 



4000m - Dmi. 



(5) 



where z is the model depth. Once sunk to the ocean floor, we assume that the 
aluminium is buried permanently. This is described by the following change of 

Alads in the bottom layer: ^ [Alads] Ig^dimentation = "[Alads] • Ws/Hsed^ where Hsed 
is the thickness of the bottom model layer. 

Adding advection and mixing, the full equations for [Alads] and [Aldiss] away 
from the boundaries are as follows: 



O"^ 

[Alads]) -g^i'^s- [Alads]) + (AV% + J/b^)[A1, 



dt 

^[Aldiss] = - k([A1:^J - [Alads]) + (AVl + UEf^)[AUissl 



adsJ , (6) 

(7) 



where the general advection (v • V) of the tracers is implicit in the full time deriva- 
tives. A and ue are the horizontal and vertical eddy diffusivity coefficients, respec- 
tively. 

The code of the NEMO model, as well as PISCES and the aluminium model, is 
free software: software that can be used, studied, and modified without restriction 
(except that published modifications must fall under the same license). It is avail- 
able at http : //www . nemo- ocean . eu/| We used version 3.1 of the NEMO model, 
specifically svn revision 1183. This can be accessed through http://forge.ipsl. 
jussieu . f r/ igcmg/ svn/modipsl/ , The aluminium specific code is available under 
the same conditions as an electronic supplement to this article: trcal_sed.F90 de- 
scribes aluminium deposition into the ocean and into the sediment, and trcal_rein . F90 
describes the internal processes of ad- and desorption. 



267 2.1.3 Input fields 

268 The dust deposition field was taken from the output of the INCA model, an at- 



269 mospheric dust model (Whitehead et al. (1998)). INCA is the aerosol module of 

270 the LMDzT atmospheric model ( Schulz et al.| ( 12009 )). The resulting climatology is 
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(a) Dust de position (g m ^ yr -*-) at the ocean (b) Sediment sources of Al (arbitrary units) as 
surface jAumont et al.| ( |2008| >) implemented in PIS CES for the Fe model 

( [Aiimont and Bopp| pOOGj ) 



Figure 2: Aluminium sources used to force the model, based on the Fe input in 
PISCES, using using an Al : Fe crustal ratio of 8.1 : 3.5. 



Basin: 


Absolute (Gg/yr) 


Relative (%) 


Atlantic Ocean 


4.49 


47.6 


Indian Ocean 


2.08 


22.0 


Pacific Ocean 


1.90 


20.1 


Southern Ocean 


0.013 


0.14 


Arctic Ocean 


0.050 


0.57 


Mediterranean Sea 


0.621 


6.58 


other basins 


0.279 


2.96 


total 


9.436 


100.00 



Table 1: Amount of yearly dust deposition in each basin. The Southern Ocean is 
defined as the ocean south of 59° S, without overlap with the Pacific, Atlantic and 
Indian Oceans. 



used and described by Aumont et al. (20 08[) a nd evaluated in Textor et al. (2006). 
The surface deposition is shown in Fig. |2(a)[ Most dust is deposited in the At- 
lantic Ocean, just west of the Sahara and just west of the Kalahari Desert. Another 
important dust deposition sites is found in the northern Indian Ocean. Dust depo- 
sition east of Asia in the Pacific Ocean is much smaller, and in the South Pacific 
Ocean, the Southern Ocean and the Arctic Ocean there is almost no deposition of 
dust. The amount of dust per basin is listed in Table [l] 

According to Textor et al. (2006) the model that is used to create this dust 
deposition field is validated against high-quality observational datasets, but there 
are a number of uncertainties, because there are not enough observations available 
to sufficiently constrain any dust model globally. The spatial resolution is nominally 
2°, and the temporal resolution is one month (twelve time steps per year). In reality, 
dust deposition is highly sporadic ( Aumont et al.| ( [2QQ8 ); Rijkenberg et al. ( 2008[ )). 
We expect that dust events are highly local and there is sufficient horizontal mixing 
to dilute the Al^iss distribution, such that there is no signiffcant concentration 
change. Furthermore, our aim is to look at averages of [Aldiss] over at least several 
months, and therefore it is acceptable to use monthly averaged dust deposition 



288 ffelds. The deposition ffeld (Fig. 2(a) ) is very similar to the one used by Jickells et al. 



289 (2005), except that in our ffeld there is almost no dust deposition near Argentina 

290 and Chili. 



In Fig. 2(b) the sediment source is plotted, which is used for one of the sensitivity 
experiments. It is a function developed for Fe ffuxes and depends on the sea ffoor 
topography (which accounts for the degree of oxygenation of the sediments), as 
described by [Aumont and Bopp (2006). The Al : Fe stoichiometry used is 8.1 : 3.5, 
the same as for the dust deposition ffeld. 
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Experiment: 


surf. diss. (%) 


subsurf.diss. 


kd (10^ 1/kg) 


K (yr ^) 


sediments 


Gehlen et al. (2003) 


3 


no 


1 


10000 


no 


Han et al. (2008) 


5 


no 






no 


Reference experiment 


5 


no 


4 


10000 


no 


Doubled dissolution 


10 


no 


4 


10000 


no 


With subsurface diss. 


5 


yes 


4 


10000 


no 


With ocean margins 


5 


no 


4 


10000 


yes 


Halved part, coeff. kd 


5 


no 


2 


10000 


no 


Slow equilibration 


5 


no 


4 


100 


no 



Table 2: Overview of reference and sensitivity experiments. Deviations from the 
reference experiment are in bold. 



296 2.1.4 Simulations 

297 Our goal in this study is to better understand how the assumptions made in our 

298 model impact on the distribution of Aldiss- Specifically, we aim to examine the 

299 importance of different Al cycle processes and external Aldiss sources. To that end, 

300 we perform five sensitivity experiments (see Table |2| with parameters different from 

301 the reference experiment. 

302 For the reference experiment we use a five percent dissolution of the Al fraction, 

303 which is within the known range and close to the values used by [Gehlen et al. 



( |2003[ ) (3%) and |Han et~al] ( |2008| ) (5%). For the partition coefficient kd we use 



4 • 10^ 1/kg. This quantity signifies the amount of Al that can adsorb onto biogenic 
silica particulates, which means that if kd is to be increased, Aldiss is exported 
in larger quantities by scavenging. Our value is four times higher than that of 
Gehlen et al.| ([2003), but much smaller than the one used by K. Ori ans and J. 



309 JMcAlister ( 5 .6 • 10 ^ 1/kg, used for shelves as well as open ocean regions, McAlister 

310 [and Orians| ( |201l[ ))p] The first order rate constant hc^ which signifies how quickly 

311 the equilibrium between Aldiss and Alads is established, is not constrained by the 

312 literature. All these parameters for the reference experiment are chosen in such a 

313 way that the model simulates the Aldiss distribution from the observations of the 



selected cruises (see Section 2.2) reasonably well, while all the parameters remained 
within literature ranges. 

After 600 yr the reference experiment is forked into several sensitivity simula- 
tions. These sensitivity experiments are run, along with the reference experiment, 
for another 1300 yr or more to a steady state. The total model time for the reference 
simulation is 2200 yr. An overview of all key parameters for the simulations is given 
in Table [2j For all experiments we use the OPA physical fields and atmospheric 



321 dust deposition as described in Section 2.1.3 In all simulations, except for the third 



which tests the sensitivity of water column dissolution, Al is dissolved only in the 
surface ocean. In the third experiment Al is dissolved in the water column as well 
with the fraction 10~^ • e~^/^° of the dust deposition field, where z is the depth in 
km and zq = 1 km. The extra term added to Eq. [rjis a • e~^^ ^^Va\/ Hk^ where Hk 
is the thickness of layer k. 

The sensitivity simulations can be divided into two types of sensitivity experi- 
ments: the first three test the sensitivity to Al sources (amount of dust dissolution 
into the surface ocean and into the ocean water column, and the inclusion of ocean 
sediment dissolution) and the last two test the sensitivity to the internal Al cycling 
(lower partition coefficient kd and lower first order rate constant n). 



-•^We did not use the kd value from [ Gehlen et aL] ( |2003 ]>, because initial experiments showed 
that this would leave too much Aldiss in the ocean surface, especially at high latitudes. 
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Cruise 


Research vessel 


Year 


ARK XXII/2 


Polarstern 


2007 


ANT XXI V/3 


Polarstern 


2008 


64 PE 267 


Pelagia 


2007 


64 PE 319 


Pelagia 


2010 


64 PE 321 


Pelagia 


2010 


JC057 


James Clark Ross 


2011 



Source 



Middag et al.'("2009 
Middag et al. ( ^201 f[ [2012 ) 
Middag et al. 
Middag et al. 
Middag et al. 
Middag et al. 



Nr. of data 



Arctic 

Southern Ocean 
North-east Atlantic 
North-west Atlantic 
Centre-west Atlantic 
South-west Atlantic 



Subtotal used primarily for detailed comparison and optimisation of the model 



Compilation of selected other observations for global ocean comparison: 



IOC96 


Knorr 


1996 


Central-south Atlantic 


M 60 


Meteor 


1982 


North-east Atlantic 


IRONAGES III 


Pelagia 


2002 


North-east Atlantic 


EUCFe 


Kilo Moana 


2006 


Equatorial Pacific 


MC-80 


Thompson 


1980 


Pacific 


VERTEX-4 


Wecoma 


1983 


North Pacific 


VERTEX-5 


Thompson 


1984 


North Pacific 


KH-98-3 


Hakuho-Maru 


1996 


East Indian 



Vink and Measu res (|2001|) 
Kremlmg|"|T9M 



Kramer et aT 



We 



et al 



:004 



Urians and BrulanS (19861) 
(1986) 
( T986| ) 



Orians and Bruland 



Orians and Brulan g" 
Obata et al. ] (|2007| 



Subtotal compilation of selected other observations for global comparison 
Grand total of all dissolved Al values used in this study 



Table 3: Observational data used for comparison with model. For positions see Fig. [3j 



2.2 Observations 

The recent IPY-Geotraces observations in the Arctic (Middag et al. (2009)), North- 



east Atlantic (Middag (2010), Chapter 5), West Atlantic (R. Middag, in prepara- 
tion) Oceans and the Atlantic section of the Antarctic Ocean ( [Middag et al.| ( [2QTT )) 



are used for a detailed comparison and optimisation of the model parameters. See 
the upper part of Table [3] for these datasets and Fig. [3] for the station positions. 
These datasets comprise overall 3455 individual data values for dissolved Al. All 
values have been verified versus international reference samples and their consensus 
values of the SAFe and Geotraces programmes, see supplementary material S-1. 
Moreover, during cruise 64 PE 321 (2010) excellent agreement was obtained for a 
complete vertical profile of dissolved Al at the BATS station between a previous 
occupation and dissolved Al data in 2008 by the US Geotraces group, see supple- 
mentary material S-2. This large amount of observations give the possibility to 
compare the model with deep ocean observations of [Aldiss], allowing us to validate 
the model in the deep ocean and to study the global Al cycle in more detail. 

For a worldwide global ocean comparison one has to rely on data that was col- 
lected in the era before the reference samples of SAFe and Geotraces were available. 
Inevitably the definition of criteria for selecting such previously published datasets 
is less rigorous, see Electronic supplement S-3 for the criteria used for each of the se- 
lected datasets. The selected datasets are listed in Table[3j lower part, and positions 
are shown in Fig. |3] 

3 Results 

3.1 Reference simulation 

Our reference experiment is performed for a spin-up period of 2200 yr. The resulting 
total ocean Al budget in the reference simulation is around 7 Tmol (1 Tmol = 10^^ 
mol) . Already after 600 yr the total Al distribution is more or less in a steady state 
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100°W 0° 100°E 



Figure 3: Stations of the cruises: M 60 [yellow]; MC-80 and VERTEX [red]; IOC96 
[green]; IRON AGES III [blue]; KH-98-3 [pink]; EUCFe [black]; Geotraces-NL [pur- 
ple]. See Table [3] for an overview with references and the number of observations. 



as shown by Fig. [4j where the total integrated aluminium (dissolved and adsorbed) 
is plott ed a gainst time. Therefore, from 600 yr our sensitivity experiments (see 
Section 3.2) are forkedj^ 

Fig. 5] shows the Aldiss surface concentration of the reference experiment. As in 
all subsequent plots, yearly averages are shown. There is a seasonal cycle in [Aldiss] 
in our model, but it is small relative to the spatial trends. The observed Aldiss 
concentrations from our Al data compilation (see Section 2.2) are plotted as circles 
over the model Aldiss concentration in this and later figures. A comparison of the 
model with the observations will be done in Sect ion [3. 1.1 1 Consistent with the dust 
deposition field (Fig. 2(a) ), the largest modelled Aldiss surface concentrations are in 
the central Atlantic Ocean, with values between 30 and 40 nM (1 nM = 10~^ mol/1) 
near 20° N, decreasing northward to values in the order of several nM north of 60° N. 
The concentration of Aldiss between 40° N and 60° N is larger than expected based 
on the dust deposition field (Fig. |2(a)| ). This is because the Gulf Stream and the 
North Atlantic Current transport Aldiss northward before most of it is scavenged. 
As we will see in Section |4.2.1| and 4^, most Aldiss is scavenged before it passes 
Iceland. The result is that [Aldiss] iii the surface waters of the Arctic Ocean is in 
the order of 1 nM. On the other side of the globe, south of 40° S, the Al concentration 
is even lower, less than 1 nM. This is to be expected on the basis of [bSi02] (Fig. 
[6| and the dust deposition fiux (Fig. |2(a)"] ). In the Pacific Ocean the concentration 
is relatively low as well, except near 20-35° N where a wide band of high [Aldiss] 
can be found. This maximum arises from the large dust deposition along the Asian 
continental margin between 20° N and 60° N. North of 40° N, the concentration of 
bSi02 is high, as can be seen in Fig. [6j Therefore, in that region the dissolved Al 
is quickly depleted by the high bSi02 concentration in the surface, while south of 
40° N the North Pacific Current transports it eastward into the central North Pacific 
Ocean. 

In Fig. [7| the modelled Aldiss concentration of the reference experiment is shown 
at selected depths (100 m, 500 m, 1000 m and 3000 m) and in Fig. [s] the Geotraces 
cruise section in the West Atlantic Ocean. This section was calculated from the 



■^The relevant tracers of the raw model output can be found at |http://data.zkonet .nl/index"7| 
|php?page=Project_view&id=29l6&tab=DatcLsets[ 
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^•0 1 \ \ 1 1 1 1 1 1 1 

200 600 1000 1400 1800 2200 

Time (yr) 

Figure 4: Total Al budget ([Alads] + [Aldiss]) (Tmol) in the world ocean, plotted 
against the model years 




Figure 5: Dissolved aluminium concentration (nM) from the reference experiment 
at the surface ocean, with respective observations plotted on top of it. Note that 
the scale is expanded for low concentrations. 
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100°W 0° 100°E 



Figure 6: Biogenic Si concentration (nM) in the PISCES model, at the surface, 
yearly average. 




1000 m 3000 m 

Figure 7: Dissolved aluminium concentration (nM) from the reference experiment 
at several depths, with respective observations plotted on top of it. 



14 



Depth (m) 



1000 ' 



5. o o 



3000 



8 o o 8 

D O O OOOO 



00 OOOO 000 

o coo OOOO OOOO 000 



O CDO O 



8 8 8 E 

o c 

OOOO 



DOOOOOOqO O O O O 00 00 00000000000 COO OOOO OOOO 000 OOOOOO OOOC 

00000000 O O O O 00 00000000000 CDO OOOO OOOO 00 OOOOOO OOOC 

2000 —OOOOOO o o o o o o o o o o 00 00000000000 CDO OOOO OOOO 000 oooOoo 000 



4000 -0 O O^ 
— o 00 o^^' 

5000 ^ " " T ^ * I 



'arii 



■000 CDO 0000 OOOO 000 OOOOQ- 



0000 coo O O c 



40°S 



20°S 



20°N 



^" 1 

40°N 



60« 



^ldfss(nM) 



Figure 8: Dissolved aluminium (nM) from reference experiment along the West 
Atlantic Geotraces cruise track. Observations are plotted on top. 



388 three-dimensional model data by converting the 0RCA2 gridded model data to a 

389 rectilinear mapping and after this interpolating the rectilinear data onto the cruise 

390 track coordinates. In the Atlantic Ocean there is a maximum of [Al^iss] from about 

391 30° S at 2 km depth, getting higher northward (Fig. |8|. According to the model, 

392 deep Aldiss concentrations are low in the Southern Ocean, the Arctic Ocean and the 

393 Pacific Ocean (Fig. [t]). 



394 3.1.1 Comparison with observations 

395 In Figs, [sj [t] and [s] observed Aldiss concentrations (see Section |2.2[ ) are plotted as 

396 circles over the modelled results. The precision of the measurements performed 

397 during the Geotraces cruises at a 1 nM level and lower (polar oceans) was about 4% 

398 and at a 20 nM level (West Atlantic Ocean) about 2% (|Middag| ([2QTq1) , p. 25). These 



399 observations have a standard deviation that is smaller than can be distinguised 

400 based on the colour bar. Older observations generally have lesser precision, but of 



401 sufficient quality (see Section 2.2). Therefore the model results can in principle be 

402 compared with the observations without regard of measurement errors. 

403 Fig. |5] shows that the model and the observations reveal similar patterns in 

404 the surface ocean: high concentrations in the central Atlantic Ocean and low con- 

405 centrations in the polar oceans and the Pacific Ocean. However, details in the 

406 modelled Aldiss distribution do not agree with the observations. For instance, the 

407 large concentration in the Pacific Ocean between 20 and 35 °N is not visible in the 



408 observations of Orians and Bruland (1986). Also between 45 and 65°N in the West 



409 Atlantic Ocean the concentration of Aldiss in the model is overestimated compared 

410 with the Geotraces observations near Greenland (R. Middag et al., in preparation). 

411 This suggests that Al is not sufficiently scavenged in this area. This is further dis- 



412 cussed in Section 4.2 The model underestimates [Aldiss] near Brazil. This can be 

413 related to errors in ocean currents (Fig. [T]), as the transport of low [Aldiss] water 

414 might be too high in the Brazil current. In reality this current is narrower than 

415 the 2° resolution of the model. Another possibility is that the dust deposition ffeld 

416 is not realistic in this area. Furthermore, the observations generally have more 

417 spatial variability than the modelled [Aldiss]- Many of these differences depend 

418 on several model boundary conditions (dust deposition distribution, velocity field, 

419 particle sinking speed) as well as temporal variability, because the model results 

420 plotted are yearly averages, while observations are done throughout the year in 
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Figure 9: The quantity in the coloured area is the model [Aldiss] (nM) at the 
Geotraces West Atlantic cruise section. The contours represent the model Atlantic 
Overturning Stream Function (Sv). 
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precise locations. 

Fig. [8] shows that similar patterns are present in the West Atlantic Ocean in 
the model and the observations. There is a clear pattern of North Atlantic Deep 
Water (NADW), which is mainly visible between 40° S and 20° N in both the model 
and the observations. In Fig. |9] the Atlantic Overturning Stream Function (OSF) 
is plotted over the model results, showing that the patterns of the OSF and [Aldiss] 
coincide. In the Southern Ocean, low concentrations of Aldiss in the upper km are 
clearly penetrating from the Southern Ocean northward to at least 20° S around 
500 m depth. This is a mixture of Antarctic Intermediate Water (AAIW) and 
Subantarctic Mode Water (SAAMW). From the Southern Ocean near the bottom, 
there is also Antarctic Bottom Water (AABW) flowing northward, visible in the 
observations and the model. 

The similarity between the model and the observations lessens in the deeper 
North Atlantic Ocean, where according to the observations [Aldiss] increases with 
depth (for depths below 800 m), while in the model there is a decrease with depth 
below 1.5 km. Besides this general pattern of increase of [Aldiss] with depth in the 
observations, a very high concentration of Aldiss is present between 45 and 50° N 
near the sediment, which enhances the dissimilarity between the model and the 
observations. This problem will be discussed further in Section [43) 

Furthermore, in the southern hemisphere in the model the maximum of [Aldiss] 
is at about 2 km depth, indicating a southward flow around 2 km depth. In the 
observations this maximum is at 3 km depth, indicating a southward flow at 3 km 
depth. Indeed, according to the hydrographical analysis during the West Atlantic 
Geotraces cruise (Van Aken et al.) the flow is at 3 km depth. As can be seen in Fig. 
[9] the southward flow is too high in the water column compared to the observations. 

Another feature which is not captured by the model is the aluminium profile in 



the eastern Arctic Ocean, which is shown in Fig. 10 It shows that the model [Aldiss] 
increases with depth from a minimum near 150 m to a maximum near the bottom 



by a factor of 3 (and is in the order of 2 nM, see Fig. 10). The observations show 
an increase as well (as also presented in Fig. [t]), but here it is far more pronounced 
with a small minimum at the surface (less than 1 nM) to a large maximum near the 
abyss (around 18 nM) ( [Middag et al.| (2009)). Thus in the Arctic Ocean the vertical 
change in [Aldiss] in the model is by far not as pronounced as in the observations. 
The absence of sediment sources of Al might play a role in this. The results of a 
simulation with sediment sources will be discussed in Section [X2I 
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Figure 10: [Aldiss] (nM) profile in the eastern Arctic of reference experiment (blue 
dashed line) and the observations of IPY-Geotraces-NL observations (black solid 
line) at 213.6°E, 87.03°N 
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Section |2.1.4| described how our sensitivity simulations are set up. The reasons 
why we performed these specific experiments are closely related to the problems we 
encountered while interpreting the results from the reference experiment (Section 
3.1|). After a brief summary of two of these problems and associated sensitivity 



simulations, we will present the results of these simulations]^ 

How can it be that in the Atlantic Ocean from 800 m depth to the bottom of the 
ocean there is a decrease in the Aldiss concentration in the model, while according 
to the observations there is an increase (see Fig. [8|? The only transformation 
mechanism in our model is reversible scavenging, which should give a profile where 
[Aldiss] increases with depth, and the only source is dust. However, the model 
biogenic silica concentration, its sinking speed, the physics and the first order rate 
constant K: are all relevant for the vertical distribution of Aldiss as well. Since 
the biogenic silica concentration and the physical fields can both be assumed to 



be reasonably realistic (see Section 2.1), no sensitivity simulations concerning these 



fields are described in this paper. We do perform one simulation with a different first 



order rate constant (Section 3.2.5) to investigate where in the ocean this parameter 
infiuences the Aldiss distribution. The value of k is significantly decreased in this 
simulation (this parameter is very poorly constrained by observations). 

Concerning the strongly depth-increasing profile of Aldiss as observed in the 
eastern Arctic Ocean (Fig. 10), we perform a simulation where sediment ocean 
jnargi ns are added as a source of Aldiss, following the approach of |Aumont and 
[Bopp (2006) and assuming an Al : Fe fraction of 8.1 : 3.5 in the sediments (thus we 
assume sediment Aldiss supply is coupled to Fe). The result of this simulation is 



^The relevant tracers of the raw model output can be found at |http : //data . zkonet . nl/ index . | 
|php?page=Project_view& id=2916&tab=Datasets 
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(a) World ocean surface difference (%) 



(b) West Atlantic section difference (%) 



Figure 11: The relative difference of the dissolved aluminium concentration between 
the experiment with an increased dissolution (10%) and the reference experiment 
(5% aluminium dissolution). The scale is in percentages, of which the largest part 
from 95% to 100%. 



discussed in Section 3.2.3 Furthermore, an experiment is performed where we have 
halved the partition coefficient from 4 • 10^ 1/kg to 2 • 10^ 1/kg, which is a way to 
evaluate the importance of relative partitioning between the dissolved and adsorbed 



483 fraction. This experiment will be discussed in Section 3.2.4 
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3.2.1 Increased dust dissolution fraction 

To examine the impact of greater dust dissolution on the cycling and distribution of 
Aldiss, we increase the surface dissolution percentage of Al from 5 to 10% ('doubled 
the resulting change in [Al^iss] in the surface ocean as well 



487 dissolution'). In Fig. 



11 



as in the West Atlantic Ocean (along the Geotraces cruise track) is presented. The 
[Aldiss] is increased with a factor of two in most locations in the ocean. 

The doubling can be explained when the model equations (Eqs.|2]and[3| are con- 
sidered. Consider a one-box model in which equilibration is instantaneous {tz = oo). 
If we increase dust dissolution everywhere with a factor of two and wait until steady 
state, we also get a doubled particulate Al burial. This means that [Alads] must 
be doubled in the ocean since the sinking/burial speed does not change. Since in 
steady state [Alads] is proportional to [Aldiss] (Eq. [3|, [Aldiss] must be doubled as 
well. Hence a doubled [Aldiss] is expected. In other words, the total Al budget is 
linear with the input. (See for instance [Broecker and Peng| ( |1982 ) for background.) 
This means that the amount of dissolution in the surface ocean has an effect ev- 
erywhere in the ocean, i.e. the effect is global and not restricted to dust deposition 
sites. 



The Al budget is not doubled everywhere. In Fig. 11(b) we can see that in the 
coastal upwelling region near Chili the increase is only around 95%. This is because 
the model is not completely spun up, as can be seen in Fig.[l2j so that the increased 
dissolution of Al from dust has not reached the deep Pacific Ocean yet. Therefore 



the percentages in Fig. 11 indicate how close the subsequent sensitivity simulations 
are to steady state in different regions in the ocean. 
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3.2.2 Water column dust dissolution 

For the reference experiment we included dust dissolution only in the surface ocean, 
to assume that Al from dust only dissolves in the surface ocean. During sinking of 
dust more Al might be dissolved. To test the effect of water column dissolution we 
performed a sensitivity simulation with water column dissolution included. 

Gehlen et al. (2003) have included Aldiss input at the surface only, but they do 



compare different dust deposition distributions, as do Han et al. (2008). The latter 



also dissolve dust below the surface layer, but no sensitivity experiments have been 
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Figure 12: Total Al budget (Tmol) in world ocean during spin- up of 1490 years 
for increased dissolution sensitivity simulation (solid black line; forked from the 
reference simulation at year 600). The dashed blue line signifies the Al budget of 
the reference simulation. 




(a) World ocean surface difference (nM) (b) West Atlantic section difference (nM) 



Figure 13: Difference in [Aldiss] (nM) between the simulation with water column 
dissolution with the reference simulation (without water column dissolution) 



presented for water column dissolution. None of our sensitivity experiments, as part 
of an OGCM, have been published before. 

Our simulation with water column dissolution resulted in slightly higher values 
of [Aldiss] of up to 0.7 nM more in some places. Fig. 13 compares [Aldiss] in the 
water column dissolution experiment with the reference experiment (Figs, [sj |8|. 

On first sight, north of the equator, water column dissolution moves the modelled 
[Aldiss] away from the observations. This is clear if we consider the decrease in 
[Aldiss] with depth in the reference simulation, while in the observations [Aldiss] 
increases with depth. The simulation with water column dissolution makes the 
modelled decrease even more pronounced. Also the high [Aldiss] in the upper 500 
m is only slightly improved, and not at precisely the right location: in the upper 
500 m the model especially underestimates [Aldiss] between 20 and 35°N. Therefore, 
water column dissolution does not improve the model results. 
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Figure 14: Experiment with sediment input: [Aldiss] (a,b) and its difference with 
the reference run (c,d) 
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It is possible that dissolution of dust is not the only way in which aluminium 
enters the ocean. We excluded river input and hydrothermal vents as a significant 
sources of Al (see Section [T]). But one other source that might be significant, and 
might also help our model to produce more realistic [Aldiss] (e.g. in the Arctic 
Ocean), is redissolution of Al from sediment resuspension. Therefore, in one of 
the simulations, we have included a rough approximation of Al margin sources, 
as described in Section 2.1.3 The resulting [Aldiss] and its difference with the 



reference experiment is plotted in Fig. 14 We can clearly see the infiuence of shelf 
sediment sources on the aluminium distribution. Especially in the Arctic Ocean, 
the Indonesian Archipelago and near the east coast of North America, dissolved 
aluminium is increased substantially when sediment input is added to the model. 
This is to be expected, since the largest sediment source is on the shelf areas (see 
Fig. |2(b) for the added Al source). 

Let's now look at the West Atlantic cruise track. Fig. 14(c)| represents the 
difference between the sediment source experiment and the reference simulations 
at the ocean surface. There is a clear increase at and near the surface between 20 
and 40° N in the West Atlantic Ocean. The cruise section difference plot in Fig. 
|14(d)| shows a large increase in the NADW. So it seems that high [Aldiss] can be 
explained by not only dust deposition but also sediment sources. This is consistent 
with the explanation of Middag et al. (in preparation), even though their arguments 
mostly concern deep sediment resuspension and to a lesser extend margin sediments. 
The sediment source in our model is mostly from near-shore sediment resuspension 



(Fig. 2(b)), in line with the findings of Mackin and Aller (1986). Of course, a deep 
sediment source might also explain the high [Aldiss] near the sediment at 45-50° N, 
and it may therefore be important to include such a source in a future model study, 
since our current results suggest that margin sediments do not effectively reproduce 
this maximum. 
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Figure 15: [Aldiss] (nM) profile from a station in the eastern Arctic (ARKXXII/2): 
observations (black solid line), the reference simulation (blue dashed line), and 
sediment sensitivity simulation (red dash-dotted line), all at 213. 6°E, 87.03°N 



556 The concentration in the model gets significantly higher in the Arctic Ocean 

557 which is bad for depths between the surface and a few km, but better for deeper 

558 waters, where a much higher concentration of aluminium was observed. This is 



559 shown in the profile plots of Fig. 15 Sedimentary input has as an effect that [Aldiss] 

560 increases everywhere in the Arctic Ocean and especially with depth, giving a profile 

561 that has the same form as the observations, suggesting the potential importance of 

562 margin Aldiss sources herein. However, the modelled [Aldiss] is too large compared 

563 with the observations. Potential reasons for this are described in Section 1431 

564 3.2.4 Decreased partition coefficient 

565 One of the key parameters in our model is the partition coefficient kd. It regulates 

566 the amount of Al that can be adsorbed by bSi02 . To see what effect this parameter 

567 has in our model setup, we performed an experiment in which we decreased kd by 

568 a factor of two. 

569 Fig. [16] shows that a decrease of kd leads to the most significant relative increase 

570 of [Al 

diss] ill the Southern Ocean. The reason for the large increase in [Aldiss] iii the 

571 Southern Ocean, and also a significant amount in the northern seas, is that these 

572 are locations of high diatom production, which results in large [bSi02]. As can be 

573 seen from Eq. [3) the equilibrium concentration of Alads is proportional to kd and 

574 [bSi02], hence kd is important in regions with high [bSi02]. This will be further 

575 explained by looking at the Al budget and timescales in Section [44] 



577 
578 



3.2.5 Decreased first order rate constant 

The second most important scavenging parameter is the first order rate constant K:. 
This parameter should have a dynamic effect on the distribution of Aldiss, since it 
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Figure 16: Experiment with a decreased kd with (a,b) [Aldiss] (^M), and (c,d) the 
relative difference with the reference experiment. 



describes how quickly Aldiss and Alads equilibrate (see Eq. [3|. To test this in our 
model, we decreased from 10^ to 10^ yr~^. 

As can be seen in Fig. 17, surface [Aldiss] increases significantly compared with 
the reference experiment (compare Fig. [s] with Fig. |17(a)] ). Fig. 17(d) shows the 
relative change in [Aldiss] in the West Atlantic Geotraces section. From the equator 
northward below 2 km depth an increase of aluminium is visible. There is also 
an increase near the deep sediment between 45 to 50° N, but it is too small to 
explain the elevated concentration in this area found in the West Atlantic Geotraces 
observations (see bottom right in Fig. l8|. 



4 Discussion 



4.1 Comparison with Gehlen et al. (2003) 



590 Since our model is very similar to the model of Gehlen et al. (2003), we expect a 



good agreement between our results and theirs. Even though several features in the 
two model results are the same, there are also noteworthy differences. These can 
arise because of differences in the ocean physics, the dust deposition field, as well 
as the biogeochemical model. 

For the velocities the output of the NEMO model on the 0RCA2 configuration 



is used in our model. On the other hand, in the model of [Gehlen et al. (2003) 



a climatology from the Hamburg Large-Scale Geostrophic (LSG) OGCM is used 
( Maier-Reimer et al. ( 1993| )). The basic characteristics of the physics of the two 
models are the same, but a few features, like the Mixed Layer Depth (MLD) and 
the Atlantic MOC, are significantly different. Even though LSG is a low resolution 
model, the depth of the meridional southward return fiow of the Atlantic MOC is 
almost 3 km deep, which is realistic. On the other hand, in the OPA model (with 
the 0RCA2 configuration) this return fiow is much more shallow (a bit over 2 km 
deep), as shown in Fig. [9] A shallow return fiow is typical for low resolution models 
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(a) World ocean surface concentration (nM) (b) West Atlantic section concentration (nM) 



Depth (m) 




(c) World ocean surface difference (%) (d) West Atlantic section difference (%) 



Figure 17: [Al^iss] (nM) in the slow equilibration experiment {hi = 100 yr ^) and 
the relative difference with the reference experiment = 10^ yi"~^). 



(e.g. [Dutay et ah] ( |2002D ). Furthermore, the AABW in LSG goes to 15°N, while 
in our model it goes all the way to 40° N. Overall, the physics of the two models 
are comparable, but we should keep in mind that the overturning in the 0RCA2 
configuration of OPA is too shallow and AABW goes too far north, compared with 
the LSG model. 



In both the work of Gehlen et al. (2003) and in our work, the highest concen- 



tration of Aldiss is found in the surface model layer (Fig. [5| in the Atlantic Ocean 
around 20° N. In both models the Aldiss concentration is low in the Southern Ocean 
and in the Pacific Ocean. However, there are distinct differences of [Aldiss] in the 
surface ocean in the models. Several features visible in the observations are better 
captured by the model of Gehlen et al.| ([2003 ) , while other features are simulated 
better by our model. These difference are likely to be due to the different dust 
deposition fields used in the model of [Gehlen et al. (2003) compared to the one we 
used. 



et al 



In the eastern Atlantic Ocean around 20° N near the coast of Africa, Gehlen 
Aldiss concentration of over 



(2003) simulate 



300 nM, or around 50 nM 
The concentrations in our model and the 



for a different dust deposition field, 
observations, however, are generally around 30 nM. 

In and near the Indonesian Archipelago [Aldiss] ranges from 1 to 20 nM in both 
models, while the observations of Obata et al."(2007) and'Slemons et al. (2010) are 
much more homogeneous, about 10 nM. This might be because the ocean currents 
in both models are not realistic because of a too low resolution of our model for 
this region. Possibly other effects in the real ocean play a role, like sedimentary 
input from the Indonesian Archipelago and east of Indonesia (see Fig. [2(b)[ for the 
relative importance). 

The vertical meridional section in Fig. [8] through the West Atlantic Ocean shows 
a reasonable correspondence between the two models. The concentration at the 
surface is high in both models, it roughly decreases with depth in the Atlantic 
Ocean. In the Southern Ocean aluminium concentrations are low in the whole 
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column, but especially in the surface (Fig. [sj very left). However, in our model 
between 10° S and 45° N there is a significant decrease of [Al^iss] with depth, which 
is not visible in the model results of Gehlen et al. ( 2003) where below 2 km north 
of the equator [Aldiss] is rather homogeneous in both dimensions along the cross- 
section. The model of Gehlen et al.| ( |2003 ) shows North Atlantic Deep Water 
(NADW) in the Aldiss distribution, spreading southward at a depth of almost 3 km, 
as also shown by the observations (Fig.jsj but also Van Aken, in preparation). Our 
model also captures NADW, but this is at only 2 km depth, as explained before. 



4.2 Advection versus scavenging 

While Al is entering the ocean through the surface by dust deposition and leaving it 
at the bottom through burial, it is also transported by the currents. In this section 
we try to find out whether scavenging or advection is more important in different 
regions of the ocean, and how this can explain certain features in the model results. 
In our model scavenging is reversible, which we base on a priori arguments 



(Bacon and Anderson (1982); Anderson ( 2006| )) and the general increase of [Aldiss] 
with depth in observations. Therefore it is unexpected that in the Atlantic Ocean 
a depth decreasing profile of Aldiss is simulated (see e.g. Fig. |17(b')] ). One possible 
reason for this is that the bSi02 distribution is not realistic. However, the PISCES 
model with the Si cycle is used by several modelli ng studies ([Aumont and Bopp 
i|2006D; [Putay et al.| ( [2QQ9[ ); |Arsouze et al.| ( [20091 ); P^tay et al.| ( |2009D ; [Tagliabue 
et al. ( 2010|)) and has t he same spatial patterns as known bSi02 export fields (e.g. 



Sarmiento and Gruber (2006)). 



We should emphasise that we only scavenge Al by biogenic silica. Of course 
CaCOa and POC might also scavenge in reality. Looking back to our comparison 
with the observations of R. Middag et al., we found that near Greenland [Aldiss] 
was overestimated by the model. This would be an example where the addition 
of other scavengers, especially CaCOs can help in decreasing the modelled [Aldiss] 
near Greenland, and also in the western South Atlantic Ocean between 30 and 40° S. 



663 CaCOa and POC; or Lam 



662 Refer to e.g. Dittert et al. (2005); Sarmiento and Gruber (2006) for export of bSi02 



(2011) for concentrations. 



The importance of scavenging is visible from several features in our model data. 
Firstly, Fig. [t] shows that on a global scale the pattern of [Aldiss] has roughly the 
same form at different depths, which is a result of reversible scavenging, i.e., at the 
surface where [bSi02] is high, Al is mainly adsorbed onto bSi02, while at depth Al is 
mostly desorbed from bSi02. Specifically, [Aldiss] is high in the Atlantic and Indian 
Oceans north of 30° S, and much lower in other regions of the ocean. Secondly, the 
fact that the concentration in the Pacific Ocean is low at all depths, is because 
Al 

diss is removed by scavenging before it reaches the Pacific Ocean. These are just 
a few of the features in the Aldiss distribution that suggest that on a large scale 
scavenging is more important than advection. 

However, the observations and the model show that advective transport by ocean 



currents is important as well. For instance. Fig. 2(a) shows that large quantities 
of dust are deposited just west of the Sahara. However, large concentrations of 
Aldiss are not only visible directly below the dust deposition site, but even higher 
concentrations are found further to the west (Fig.[5|. Large amounts of dissolved Al 
must be advected from the dust deposition site toward Central America. Hereafter 
it seems to be advected northward until Iceland, where [Aldiss] is high despite very 
little dust deposition. Deeper in the North Atlantic Ocean the original dust signal 
decreases, and at 2 km depth in the model (or 3 km depth in the observations) 
Aldiss is transported into the southern hemisphere by the NADW. This southward 
transport can be seen clearly in Fig. [Sj but no further than 40° S. South of this 
latitude the Aldiss concentration is low. These examples indicate that also advection 
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686 plays an important role in the redistribution of Al. 

687 4.2.1 Timescales 

688 While the analysis of [Al^iss] distributions points to the combined importance of 

689 scavenging and advection in setting its overall pattern, it does not allow to quantify 

690 their relative importance. In order to assess the relevance of these different processes 

691 we define associated timescales based on our model equations. 

A priori the relative relevance of scavenging can partly be derived from the 
model equations. If Eq. |3]is substituted into Eq. [2j we get: 

^[AW.] = ^-(TOJ-[AW.]) 

= — [AldiJ-— [Alads], (8) 

692 where rtrans = l/{i^kd[hSi02]) is the typical adsorption timescale. This is the typical 

693 time necessary for Al to transform from Aldiss to Alads (or vice versa). This process 

694 depends on the amount of available bSi02 and the amount of Aldiss which can 

695 adsorb onto bSi02, giving the actual rate of conversion between Aldiss and Alads- 

696 The addition of Alads is equal to the removal of Aldiss , so rtrans is either the time 

697 necessary to add Alads, or the time to remove Aldiss- The second timescale in Eq.[8) 

698 Teq = 1/hz/is the time it takes for Alads to equilibrate to Al^^^. This process provides 

699 a stabilising feedback on the growth of [Alads] (or decrease of [Aldiss] )- 

700 Scavenging is the process of adsorption and sinking, so we need to know how 

701 fast particles sink. We use this together with the adsorption timescale to define a 

702 scavenging timescale. Since we look at the surface ocean, we might define the typical 

703 sinking speed as sinking through the first 10 m, the thickness of the upper model 

704 gridbox. However, after exporting Alads out of a gridbox by sinking, the particles 

705 sink through the mixed layer (which is generally thicker than 10 m). During sinking, 

706 the particles will partly be mixed back into the upper gridbox. Therefore to speak 

707 of a removal, it is a necessary condition for Al to sink out of (at least) the mixed 

708 layer. In our model, the sinking velocity of bSi02 and Alads is defined to be constant 

709 in the mixed layer, namely Ws = 30 m/day. Now we are easily able to define the 

710 typical sinking timescale: Tgink = ^ml/'^s, where I^ml is the mixed layer depth. 

711 The scavenging timescale can be defined as the maximum of the adsorption and the 

712 sinking timescales: 

'^scav ~ I^ax('n;rans7 '^sink) - (9) 

713 This signifies how fast Al is exported from the mixed layer. Except for a few places 

714 where the mixed layer is very deep, like the centre of the Labrador Sea, the time it 

715 takes to sink out of the mixed layer is generally very small compared to the time it 

716 takes the Al to adsorb onto a particle (rsink ^ ^trans) and thus almost everywhere 

'^1'' '^scav ~ 'Hirans- 

718 A typical advection timescale is defined as follows: 

Tadv = L/V . (10) 

719 Here L is the typical length scale, defined as the horizontal diameter of a gridbox, 

720 and V is the horizontal speed. Since the vertical velocity component is very small 

721 (W <C {Dml/L)U), only the horizontal components need to be considered. The 

722 meaning of this timescale is that within a time Tadv a dissolved Al particle (or any 

723 other non-buoyant tracer) is advected out of a gridbox (if it is not scavenged before). 

724 The residence time can be defined as the minimum of the scavenging and the 

725 advection timescales. It is the typical time that a Aldiss particle stays within a 

726 volume box of horizontal resolution (about 2° x 2°) times I^ml- This quantity 
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Figure 18: Residence time in months of Al in a box of horizontal resolution times 
I^ML, or Tres = min(radv7 '?"scav) • It represents the average time an Al^iss particles 
spends in this box, considering only scavenging and horizontal advection. 



727 is presented in Fig. [Tsj The oligotrophic gyres in the Atlantic, Pacific and Indian 

728 Oceans are clearly visible. In the centre of these gyres it takes much longer than one 

729 year before Al is exported out of a volume box, either by advection or by scavenging. 
The modelled Aldiss distribution of Fig. [5] shows in the Atlantic Ocean the largest 
values between 10 and 30° N. This can be partly explained by the large dust input 

732 between 5 and 20° N, where Aldiss is removed relatively fast by advection but dust 

733 input keeps [Aldiss] large, and partly by the large residence time between 20 and 

734 30°N, where there is no large dust input but Aldiss simply stays there for a long 
time, since there is no bSi02 present for scavenging and advection is very small. A 
similar argument can be given for the high [Aldiss] near the North Pacific gyre. 

737 If we want to know which of the timescales are more important, a relative relation 

738 between Tscav and Tadv must be defined. The number for relative importance of 

739 scavenging versus advection for Al export is the following: 

T=^. (11) 



730 
731 



735 
736 



740 



748 
749 



The logarithm of this quantity is plotted in Fig. 19 In regions where T ^ 1, like 

741 north of 40° N in the Atlantic and Pacific Oceans and in the south of the Southern 

742 Ocean, scavenging is more important than advection. On large scales the large T 

743 regions coincide very well with the large bSi02 regions (Fig. [6|. Looking more in 

744 detail, we can see that T is less homogeneous than [bSi02]. For instance, in Drake 

745 Passage and in the Pacific sector of the Southern Ocean at 150°W the timescale 

746 fraction is smaller then one would expect based on the bSi02 distribution. This is 

747 because of the strong current of the Antarctic Circumpolar Current (see Fig. [T]). 
In regions where T <C 1, in the low latitudes, advection is more important than 

scavenging. One can say that if T <C 1, Al export is advection driven, while for 

750 T 1, Al export is scavenging driven. In the high T regions a one-dimensional 

751 model would be a reasonable approximation, provided that there is sufficient Al 

752 input to keep [Aldiss] in steady state. One region where this should perfectly work 

753 is the Mediterranean - there T is very large and the dust fiux is high. Several 

754 other regions like the northern seas, near Antarctica and some other coastal areas 

755 would also be good candidates for a one-dimensional model. Other regions of the 

756 ocean cannot be well described by a one-dimensional model. Specifically the Arctic 
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Figure 19: logio(T), when hi = 10^ yr~^ and kd = 4: ■ 10^ For the advection 

time Tadv the definition of Eq. 10 is used, and the scavenging time Tgcav is the 
maximum of rtrans and Tgink- 



Ocean would not be a good candidate. Even though it is a semi-closed basin, 
internal circulation has an effect of the same order as scavenging. These arguments 
must be taken into consideration when using one-dimensional models. 
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4.3 Ocean sediments source 

Including margin sediment input of Al results in overestimated values for [Aldiss] 
compared to the observations in the Arctic Ocean. Potential reasons for this will 
now be discussed. 



As described in Section 2.1.3 the Al flux is proportional to the Fe flux which 
depends on the degree of oxygenation of the sediments. This is probably not a rea- 
sonable assumption, since generally Al does not flux out of undisturbed sediments, 
but enters the water column through resuspension of the sediments (e.g. iMackin 
(1986); Van Beusekom et al. ( 1997[ )). Since there is not enough data on such resus- 
pension, we choose this simplistic parametrisation for our sensitivity experiment. 
A priori we therefore have no reason to expect an improvement in the model. The 
sensitivity experiment can only be used to get an idea of the first order effects of 
an ocean sediment source of Al, henceforth referred to as "sediment input" . 

Since the sediment input in the Arctic Ocean is treated in the same way as the 
rest of the ocean, it is surprising that [Aldiss] gets much too high (Fig.[l4|, certainly 
when considering the improvement in the West Atlantic as described below. There 
is a lack of knowledge of boundary exchange processes (e.g. Arsouze et al. ( 2QQ9[ )), 
which keeps open the possibility that the large sediment source from the margins is 
compensated by boundary exchange processes. In our model only bSi02 was used as 
a scavenger, while near the margins, other scavengers like POC and CaCOs might 
play an important role. Also scavengers not present in the PISCES model, like 



clay minerals (Walker et al. ( 1988, )) might be important for reducing [Aldiss]- The 
possibility of this boundary scavenging (Bacon et al. (1988); Arsouze et al.|(|2QQ9)) 



is convincingly shown to be present in observations (e.g. Brown et al.| (2Q1Q|)). 



Since most of the data used here is from the open ocean and not from coastal 
areas, and our model is too coarse for an analysis of the coastal aluminium con- 
centrations, it is not possible to make a good a priori estimate of the amount 
of sediment input needed for a better simulation. Since dust input and internal 
processes contain parameters which are, just like sediment input, not completely 
constrained, a change in sediment input could result in a simulation which predicts 



27 



1 




























/' 




/ 




/ 




/ 




/ 




/ 




/ 




i 











100 300 500 700 900 1100 1300 

Time (yr) 



Figure 20: Total Al budget (Tmol) in the world ocean after a partial spin-up of 600 
yr of the reference simulation. From this point the reference simulation together 
with the sensitivity experiments is run for another 1300 yr or more. The blue 
solid line is the budget of the reference simulation, the red dashed line of the slow- 
equilibration simulation and the black dash-dotted line of the simulation with a 
decreased partition coefficient. 



790 
791 
792 
793 
794 



the open ocean aluminium concentration better 
eas. Nevertheless, consistent with work by 



but does a bad job in other ar- 
p9ll); 



Moran and Moore 



Moran et al. 



(1992); Middag et al. (2011), there are clearly some areas where sediment supply of 



Aldiss appears important and a better understanding of the processes governing its 
supply is needed. 
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4.4 Internal coefficients 

For the experiments where we changed two internal parameters (kd and we can 
now analyse the simulations in a more sophisticated manner by using the above 
timescale approach. 



As can be seen from the black dash-dotted line in Fig. 20 , the effect of a halved 
partition coefficient on the whole ocean is an almost doubling of the total Al 
budget. According to Eq.js] indeed a halved kd means half as much [Al^^^]. In all 
experiments equilibration is very fast, except the one where n is decreased signifi- 
cantly. Because of the fast equilibration, with halved kd^ export will be half as fast. 
In a steady state, this would result in a doubled Aldiss content, if bSi02 would be 
distributed homogeneously. Since it is not, this does not hold for our model. At the 
dust deposition site in the Atlantic Ocean we see in Fig. [16(c)] an increase of [Aldiss] 
of less than 100% in the Atlantic Ocean near the dust deposition site. Because of 
the reduced kd more Aldiss can be transported northward. When it reaches large 
[bSi02], it is scavenged with twice the speed compared to the reference simulation. 
Since this extra Aldiss was able to reach the high [bSi02] site, the effect of the 
reduced kd is dampened and therefore the ocean Al budget is less than doubled. 

If we also want to consider advection, we should look at T, which is proportional 



to kd and [bSi02] (see Eq. 11). The halved kd results in only a small change in 
importance of scavenging relative to advection (see Fig. 19 where the logarithm of 
T is plotted). 
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A decrease in the first order rate constant hi means that equilibration goes slower 
(Eq. |2|. If dust is dissolved in the surface ocean, its conversion to Alads is slower, 
so that Aldiss is scavenged slower. Because of the resulting higher [Aldiss], more 
aluminium is transported northward by the North Atlantic Current. Fig. [6] shows 
the biogenic silica (bSi02) concentration according to the model at the surface. In 
the Atlantic Ocean, we can see large concentrations of bSi02 north of 40° N. The 
dissolved Al which arrives in this area, is scavenged by this high [bSi02]. This 
results in high concentrations at all depths around 40° N. However, because of slow 
equilibration, Aldiss can go further north before it is actually scavenged. This is 
the reason why a significant increase in [Aldiss] is visible around 60° N at all depths, 
and also further south at a depth around 2 km because this extra Al is advected 
southward via the MOC (Fig. |17(d)[ ). Also desorption goes slower, which means 
that Aldiss that is adsorbed at the surface around 40° N, will desorp slowly while 
sinking, such that compared to the reference simulation a relatively large quantity 
of Aldiss will appear in the deep Atlantic Ocean. This can be seen in the lower right 



831 of Fig. 17(d) 



According to Fig. 17(d) there is a decrease of [Aldiss] at low latitudes around 500 
m and 1 km depth. Less dissolved Al from dust in the surface ocean is scavenged 
and therefore less adsorbed aluminium is present in the deeper regions to desorb at 
depth. So in the deep ocean there are regions of increased [Aldiss] and regions of 
decreased [Aldiss]- 

If a lot of particulate aluminium sinks where there is almost no biogenic silica, 
it will desorb and the resulting Aldiss will stay there. According to the equations, 
with a very large K: this will happen instantaneously and the dissolved aluminium 
will stay there indefinitely (in the limit of zero biogenic silica and zero advection), 
i.e., Tres = oo. This really means, as discussed in Section 4.2. 1[ that in that case 
advection is much more important than scavenging. If is decreased, the particulate 
aluminium has enough time to fall out of this low biogenic silica domain before it 
desorbs. 

In Fig. 21 the logarithm of the relative importance of scavenging, T, is plotted 
for the slow equilibration experiment. It shows that everywhere advection is more 
important than scavenging. Note that the residence times are based on boxes of 
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40°S 20''S 0° 20° N 40°N 60°N 40"s 20-s o- 20-n 4o=n 

(a) (b) 

Figure 22: The quantities in the shaded area are (a) [Aldiss] (nM) and (b) [Sidiss] 
(/iM) at the Geotraces West Atlantic cruise section. The contours represent the 
Atlantic Overturning Stream Function (Sv). See Fig. [9] for the large version of (a). 



approximately 2° x 2° x I^ml- Fig. 21 can be interpreted as that with a low 



849 the aluminium distribution gets more spread out, or homogenised, throughout the 

850 ocean compared with the reference experiment. As can be seen from the thick red 
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851 line in Fig. 20 the relative increase of the Al budget is very small. The only strong 



852 increase is in the surface ocean. In the rest of the ocean, [Aldiss] is just homogenised. 

853 4.5 Al versus Si in the MOC 

854 Our modelling results in the Atlantic Ocean show that at around 2 km depth [Aldiss] 

855 decreases from north to south, while the concentration of dissolved silicon (Sidiss) 

856 increases, all the way into the Pacific Ocean, both following the Meridional Over- 



857 turning Circulation (MOC). This is shown in Fig. 22, representing the concentra- 

858 tions of Aldiss and Sidiss in the West Atlantic Ocean. On top of these figures a 

859 contour of the Atlantic Overturning Stream Function in Sv (1 Sv = 10^ m^s~^) of 

860 the physical forcing is plotted. 

861 The patterns for both [Aldiss] and [Sidiss] are in first order easily interpreted 

862 by looking at the source of the tracer and the general currents in the Atlantic 

863 Ocean. The source of Aldiss is dust dissolution in the surface of the central Atlantic 

864 Ocean. It is transported northward by the Gulf Stream and the North Atlantic 

865 Current. Before the Aldiss reaches locations where NADW is formed, it is scav- 

866 enged by biogenic silica from 40° N northward (Fig. [6| and based on our timescale 



867 analysis (Fig. 19) possibly already near 30°N. Because of remineralisation [Aldissj 

868 increases in the NADW. This is consistent with the fact that in the model Aldiss 

869 does not nicely follow the Atlantic Overturning Stream Function in the very north 

870 of the Atlantic Ocean, but rather tends to be transported into the NADW indeed 

871 near 30° N. After this the Aldiss is taken along with the NADW, reaching 40° S at 

872 around 2 km depth. To complete the picture of the modelled Aldiss distribution, 

873 low [Aldiss] SAAMW/AAIW enters just below the surface of the Atlantic Ocean 

874 from the Southern Ocean, and low [Aldiss] AABW into the Atlantic Ocean below 4 
km depth. 

The source of Sidiss is the large amount of diatoms that sink and remineralise 

877 south of 50° S, after which it is advected as AABW as clearly shown in Fig. |22(b)[ 

878 Other features visible are the SAAMW/AAIW moving below the surface of the 

879 Atlantic Ocean and low [Sidiss] NADW diluting the dissolved Si concentration in the 

880 Atlantic Ocean. The pattern of [Sidiss] is similar to the observations described by R. 

881 Middag et al., in preparation. The main difference is that the model overestimates 

882 [Sidiss] in the deep North Atlantic Ocean, and is a little bit too low in the deep 

883 South Atlantic Ocean. This is consistent with the overestimated AABW inflow in 

884 the model velocity field. In reality the lower cell of the Atlantic OSF reaches a 
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maximum of 1 or 2 Sv, while according to our physical model forcing its maximum 
is about 6 Sv. Furthermore, the AABW does not go as far north as in our model. 
Because of this the transport of any tracer following the pathway of AABW will 
be too large. Therefore [Sidiss] is overestimated in the deep North Atlantic Ocean. 



Coming back to the discussion in Section 3.1.1, because of the same reasons [Al^iss] 
is too quickly diluted in the deep North Atlantic Ocean. 

The underestimation of [Al^iss] in the deep North Atlantic Ocean is more sig- 
nificant than the overestimation of [Sidiss] in this same area. This suggests there 
might be another process playing a role. A likely candidate is an adjusted depth- 
depending sinking velocity of bSi02 and Alads which is depth-increasing in our 
current model. Recent observations (iMcDonnell et al. ( 2Q11| )) suggest that sinking 
speed is not strictly depth-increasing, but we should consider other functions with 
depth or more sophisticated aggregation methods, as for instance done in models by 



Kriest and Evans 


(1999 


); 


Gehlen et al. (2006); 


Burd and Jackson 


(2009 


), but none 


of these efforts give results consistent with the study by 


McDonnell et al. ( 


2011). 



Another possibility is that we need sediment sources more sophisticated than in 



our sensitivity experiment, as also shown by observations by Moran et al. (1992) 
Middag et al. (in preparation) and references therein. 



903 5 Conclusions and outlook 

904 The objective of this study is to come to a better understanding of the behaviour 

905 of Al, by means of simulating the Aldiss distribution in the ocean with a reversible 

906 scavenging and general circulation model, and comparing it to observations. In 

907 turn, these results can be used for further developing this model to simulate the 

908 Al distribution more precisely. A more realistic simulation could then be used for 

909 constraining dust deposition fields to more precisely know nutrient input rates. It 

910 can also be used to study further the infiuence of sediments on the water column 

911 [Aldiss] 7 e.g. by means of coupling the ocean model to a sediment model. 

912 We ran the biogeochemical model PISCES off-line, forced by a climatological 

913 velocity field with a temporal resolution of five days and a monthly dust deposition 

914 field. In one of the sensitivity experiments we included a margin sediment source, 

915 which has not been done before for Al. We are able to simulate the main features 

916 of the global Aldiss distribution in accordance with available observations. Specifi- 

917 cally we are able to simulate reasonably well the distribution in the West Atlantic 

918 compared to observations from the West Atlantic Geotraces cruises. Since our re- 

919 suits are close to observations, we can assume that the Al distribution is indeed 

920 mainly controlled by advection and reversible scavenging, with bSi02 as the main 

921 scavenger. 

922 It is possible to improve specific features of the distribution by changing certain 

923 parameters of the scavenging process, or adjusting the Al sources. Increasing dust 

924 Al dissolution results in an overall aluminium increase proportional to the increased 

925 dust factor everywhere in the ocean. Adding dissolution in the water column gives 

926 higher concentrations, especially near dust deposition sites around one or two km 

927 depth. Decreasing the partition coefficient results in a higher concentration of 

928 Aldiss everywhere in the ocean. Especially the relative increase in the Aldiss concen- 

929 tration in the Southern Ocean is large. This parameter highlights the importance 

930 of the inffuence of spatial variability in biogenic silica on [Aldiss], since kd signifies 

931 the amount of Aldiss which can adsorb onto bSi02. When we include sediments as 

932 an input to Aldiss, we see elevated concentrations especially near the margins, as 

933 expected. [Aldiss] is strongly elevated in the Arctic Ocean, which is not reconcilable 

934 with the observations, and is possibly due to the lack of oceanic-margin scavengers 

935 in the model. The model should be modified such that a realistically constrained 
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ocean margin source of Aldiss is present. 

By means of a timescale analysis based on our model equations, we have shown 
that the importance of scavenging versus advection is highly location dependent. 
The gridbox based residence time is strongly space dependent as well. Especially 
latitudinal gradients are large. This residence time varies between less than a week 
at high [bSi02] and strong advection, to many years in the oligotrophic gyres. 

Even though the most relevant features of the Aldiss distribution are captured, 
some fine-tuning of the key parameters is necessary. There is clearly a localised 
source missing in the West Atlantic between 45-50° N near the sediment. As already 
hypothesised by Moran and Moore ( |1991 ), resuspension of nepheloid layers along 
the western boundary of the North Atlantic Ocean is a source of Aldiss- A more 
detailed analysis of the locally elevated [Aldiss] in the West Atlantic near 45° N will 
be published by R. Middag et al. (in preparation). 

Another location of interest is the Arctic Ocean, where we could not simulate the 
observations. This is possibly because we miss a sediment source, or an important 
process is missing in the model, namely the biological incorporation of Al into the 
diatom's frustules fCas chetto and Wollast] ( |1979D ; [Gehlen et al.| ( |2QQ2D ). When this 
process is implemented in the model, dissolved Al could be taken up by diatoms 
in the surface ocean and released in the deep ocean when the biogenic silica is 
remineralised together with the incorporated Al, yielding a strong Al : Si correlation 
as present in the observations (Middag et al. (2009)). 

We made good progress in simulating the distribution of dissolved aluminium in 
the world ocean, and our approach confirms that dust deposition is the main source 
of aluminium and reversible scavenging is the main process in removing it. However, 
significant improvement might be possible by (1) developing a more sophisticated 
model for ocean sediment source of Al, (2) adding diatom Al incorporation, (3) 
understanding where the elevated Al levels in the North Atlantic Ocean near the 
sediment come from and incorporating these sources in the model, and (4) fine 
tuning the free parameters and constraining them by means of field studies and 
laboratory experiments to better understand the role for these processes. 
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